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Abstract 

Nonlinear integrable models with two spatial and one temporal 
variables: Kadomtsev-Petviashvili equation and two-dimensional Toda 
lattice are investigated on the subject of correct formulation for bound- 
ary problem that can be solved within the framework of the Inverse 
Scattering Problem method. It is shown that there exists a large set 
of integrable boundary problems and various curves can be chosen as 
boundary contours for them. We develop a method for obtaining ex- 
plicit solutions to integrable boundary problems and its effectiveness 
is illustrated by series of examples. 

1 Introduction 

Known examples of integrable 2+l-dimensiona|^ systems are still not nu- 
merous, but they attract considerable interest in the aspect of the theory of 
differential equations as well as of various applications. Probably the most 
noteworthy example is Kadomtsev-Petviashvili equation (KP, see [T]) 

Ut + Uxxx — QuUx = 
Wx = Uy, 

where u = u{x,Y,T), w = w{x,Y,T), 
variables, parameter a equals either 1 or i. 

^Usually dependence on two spatial and one temporal variables is implied 



1^, X, F, T are continuous 
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Another important example is presented by difference-differential model 
- two-dimensional (2D) Toda lattice (see [2j) 

vxy{X, F, n) = W{X, Y,n-l)- W{X, Y, n), (1 .2) 

where f (X, F, n) is an unknown function of two continuous variables X, Y 
and integer index n G Z, W{X, Y, n) = exp{v{X, Y, n) — f (X, Y,n + 1)). 

The phenomenon of integrability of such systems is based on existence of 
so-called commutator representations or Lax representations. Examine this 
concept on the example of KP system (11 .11) . It was shown (see papers [3], 
[1], [5]) that nonlinear equation (11 .11) constitutes consistency condition for 
two linear equations of the following type: 

which specify differential Lax operators L, A so that KP equation is written 
in the commutator form 

[dy-L,dT-A]={), (1 .4) 

where 5y = d/dY , [, ] is commutator. For this reason KP equation fll .ip may 
be referred to as integrable model which possess a large class of solutions 
important for both theory and applications. There exist soliton solutions to 
KP that can be obtained with the help of the "dressing" method (see [3], 
0)- 

The special interest from the point of view of applications is inspired 
by possibility to solve boundary problems for 2+1-dimensional models. It 
must be mentioned that not any boundary problem for KP is solvable in 
the sense of presence of large enough set of representative "nice" solutions 
like those solitons. The point is that by imposing an arbitrary boundary 
constraint one brakes the standard scheme of the Inverse Scattering Method 
(ISM) since the infinite hierarchy of integrals of motion vanishes. So we arrive 
at the question if there exist formulations of boundary problem which save 
sufficiently many integrals of motion to keep up applicability of the ISM. 
Considerable progress in this aspect was provided by papers [6], [7], [8], [9], 
|lUj . The following Definition constitutes the condition which allows to select 
boundary problems compatible with the ISM. 

Definition 1. Let nonlinear equation admit two different Lax representa- 
tions not connected via any gauge transformation: 

{^x = -4^', ^'y = L^}, {§x = = L^!}. 
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Then the boundary constraint Q{u) =0 on the boundary X = Xq is called 
integrable if there exists a linear differential operator B so that on boundary 
X = Xq function \1/ = B^^ solves the equation \Efy = L"^ for any solution \Ef 
of equation \E'y = L\Ef as the constraint fl{u) = holds. 

In papers [8], [9] the Definition 1 was, in partial, used to answer the ques- 
tion of correct formulation of boundary problems for KP equation on straight 
boundary contour (X = Xq), similar questions for 2D Toda lattice with cir- 
cular boundary contour were investigated in paper [15]. Nevertheless it was 
unclear how to determine the class of integrable (in the sense of Definition 
1) boundary problems, which curves may be chosen as appropriate boundary 
contours. 

This paper's goal is to develop a method that constructively determines 
integrable boundary problems for 2-M-dimensional models, KP equation 
(11 .ip and 2D Toda lattice (II .2p form main examples. We present an algo- 
rithm for integration of corresponding boundary problems and the procedure 
yields solutions in the explicit form. 

In section 2 we show main integrability attributes for equation (II .11) : 
Gelfand-Levitan-Marchenko (GLM) system associated with Lax representa- 
tion (II .3p and dressing method. In section 3 we rewrite equation (II .ip in 
other variables and determine appropriate boundary contour and boundary 
constraint in generic form. In section 4 the mentioned tools are used to solve 
concrete boundary problem: we write out the GLM equation and find its 
kernel that enjoys the integrability condition. Then we utilize the dressing 
method to find solutions to boundary problem in the explicit form. The last 
part of the paper applies the developed method to solving boundary problems 
for 2D Toda lattice (O]). 

2 The use of the ISM. 

Remind, following 0) 0, Zakharov-Shabat dressing method, one of ways to 
obtain explicit formulae for solutions to KP equation. Introduce operator F 
such that 

^ roc 

(F^)(x) = / F{x,z,T,Y)'i'{z)dz, (2.1) 

J —oo 

[Ot -a, F] = 0, [dy -L, F] = 0, (2 .2) 

and suppose that it admits the triangular factorization, i.e. the following 
form: 
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1 + F = {l + K+)-\l+K), (2.3) 

where 

/X C OO 

-oo J X 

(2 .4) 

Now examine operators A^: 

L^ = dl- 



X 

= -id't + u^d^ + ?,{ul + awl,, 



(2 .5) 

specified by the following constraints: 

{dT-A^){l + K) = {l + K){dT-A), 

{dy - + K) = {1 + K){dy - L). 

They satisfy the commutator relation (11 .41) . have the structure similar to 
that of operators L, A and their coefficients can be found with the help of 
equations (12 ■6p . so that functions form new solutions to KP equation via 
the following formula: 

= u + 2dxK{x,x,Y,T) (2.7) 

To find the function K(x, zj^ we need to solve GLM equation which arises 
from systems (12 ■3p and (12 .41) : 

K{x,z) + F{x,z)+ r K{x,z')F{z',z)dz' = 0. (2.8) 

J —oo 

Thus the dressing procedure consists of the following: we start from a sim- 
ple potential u{x) and solve linear equations (12 .21) on kernel F{x,z). Then 
substitute F{x,z) into GLM equation (12 .8p . find K{x,z) which via formula 
(12 .71) yields new "dressed" solution to equation (II .ip . 
Rewrite system (12 .2p in the following form: 

aFy - F,, + F,, + {u{x) - u{z))F = 
Ft + ^F^xx + 4F,,, - Qu{x)Fx - 6u{z)F, (2 .9) 

+3[a{w{z) — w{x)) — Ux{x) — Uz{z)]F = 0. 

Look for solution to system (12 .9p in the form of product F{x, z) = '${x)'${z). 
It is easy to verify that function ^(x) solves system (II .31) whereas "if^x) solves 
system 



^For the sake of brevity we omit the dependence on variables Y, T. 
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= -4§^^^ + 6m§^ + 3{u^ - aw)^, 

which we call conjugate to system fll .3p . Then substitute just obtained kernel 
F{x,z) into GLM equation fl2 .81) and deduce formula for K{x,z): 

K(x, z) = ^ ' , 2 .11 

l + I^^n^'Mz')dz'' 

which with the help of (12 .71) yields the "dressed" solution to KP. 

The separate question is how to select solutions to integrable boundary 
problems among all the dressed solutions we obtained above. This will be 
done further via Definition 1. 



3 Integrable boundary problems. 

Having the operator L— A-pair fll .31) one can rewrite it via variable transform 
Y = Y{y,t), T = T{y,t) so that integrable through Definition 1 boundary 
constraint is posed on straight contour y = yo or, in former variables, on 
curve line 

Y = Y{yo,t), T = T{yo,t). (3.1) 
In this case the system of linear equations (II .3p looks in the following way: 



= -i {^(^- - - /ii(-4^xxa= + Qu^. + 3(m, + , 

(3 .2) 

where fij are partial derivatives: 

T T 

hi = yY = 7^ 777^, Jl2 =tY - 



■J-tJ-y ~ J-t-Ly J-ylt — Iylt 

Yt Y 

/2I = UT = 777T? ^TTTT' /22 =tT = TfT^ ^7T7' ^ ^ /11/22 — /12/2I, 

J-y'^t '~ J-tiy ^tiy -'y-'t 

and conjugate system (12 .101) in new variables is obtained from (13 .21) by 
simple replacement «—)■—«: 

= -i {^(^- - + /i2(-4^.xx + Qu^. + 3(m, - aw))^ 
^t = i {^(^xx - u^) + /ii(-4§,,, + 6m§, + 3(u, - aw))^ . 
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We associate the system fl3 .21) and its conjugate one fl3 ■3p with the same 
solution u{x) of integrable boundary problem, so in virtue of Definition 1 
we demand that linear correlation \E'(x) = B(x,t)'^(x)\y=y^ holds on the 
boundary contour. Substitute this into the second equation of system (13 .Sj) . 
with the help of (13 .2p collect coefficients at functions \E', "^x, "^xx, "^xxx and 
arrive at explicit formulae for the multiplier B{x,t) : 

B{x,t) = g{t)e^ , (3.4) 
where g{t) is arbitrary function, and for the boundary constraint: 



a;A //2l\ , /21M , ^ , /; 



A(log(7(t)), + — +i^ + Qah^w- 



3 

21 



6a l/ii/ a lOSa^/fi 



0. (3 .5) 



y=yo 



We must also demand that the following conditions on the variable transfor- 
mation hold: 

hi + 0, /21 ^ 0, A 7^ 0. (3 .6) 

4 Explicit formulae for solutions to integrable 
boundary problems. 

The variable transformations Y = Y{y,t), T = T{y,t) specify boundary 
contours as parametric curves (13 .ip . For example we choose the following 
simple transformation: 

Y = yt, T = |, (4.1) 

which determines the boundary contour as hyperbolic curve YT = ?/g. In this 
case the integrable boundary constraint (13 .5p is presented by formula 



X 



+ u + 6a —w 



3T Y 108 T2 



= 0, (4 .2) 

YT=yl 



for simplicity g{t) = const. 

Let the following elementary solution of the boundary problem (II .ip . 
(14 .2p be the initial step for the dressing procedure: 



^In general i? is a differential operator, here for simplicity we restrict to the case of B 
being a constant multiplier. 
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18T3 ' 18T3 ^ 36T4 648T6 ^ ^ 

Substitute solution fl4 ■3p into coefficients of linear systems fl3 ■2p . f l3 .3p and 
find solutions of these systems: 



232/0^ 



^ = Ef=i dj exp 



36T3 



48T5 



+ 



36T3 



aT 



E^i exp 

36T3 



12ar2 y 



12aT2 



48T5 



2 4 ,,2„2 

36T3 aT 



(4 .4) 



(4 .5) 



where dj, Ij, pj, qj, j = 1,2,. ..,N are constant parameters. Choose the 
simplest case = 1 which may be formally associated with one-soliton 
solution. To keep the integrability one must satisfy Definition 1, i.e. the 
following linear condition on the boundary: 

^ = B^\y=y„ S = e^. (4.6) 

So we must take solutions (14 .41) . (14 .5p of linear systems (13 .21) . (13 .3p and 
select those solutions which satisfy the condition (14 .6p . For = 1 and p = q 
this condition holds. Then substitute the appropriate functions \E' and \E' into 
the formula for solutions to GLM equation: 



K{x, x) 



1 + J^^^{z)^{z)dz 




(4 .7) 



where 



X = exp —2xp + 



a^Yylp 



1 A 



+ 8Tp^ 



3T2 18T3 

and thus in force of formula (12 .7p obtain the "dressed" solution to KP: 



(4 1 



U 



w 



2 2 



+ 



18T3 

9 9 9 

3T2 cosh^ - 



cosh^ 
V 

+ 



36r4 



(4 .9) 



648r6 ' 



where 



r = —px + 



4T2 



2 4 



1 



36T3 +4Tp=^ + ^M2p). 
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Note that putting yo = we obtain straight boundary contour Y = 0, hnear 
phase T and the solution f l4 turns into an ordinary sohton solution of KP 
and Korteweg-de Vries equations. 



5 Boundary problems for discrete problems. 

The method we developed above can be successfully applied to solving bound- 
ary problems for integrable 2+1-dimensional chains. Examine this procedure 
on the example of 2D Toda lattice (11 .2p . For convenience rewrite it in other 
notations: 



where w{X, Y, n) = exp{u{X, Y, n) — Y,n + 1)). 

This model was proposed in paper [2] together with appropriate Lax 
representation and integrals of motion. The nonlinear chain (15 .11) presents 
the consistency condition for the following pair of linear equations: 



where are difference Lax operators. Authors of papers [TT], [12] used 

this fact to develop approach to systems like (15 .11) from the Inverse Scatter- 
ing Problem viewpoint. Papers [13], [H] presented the dressing method for 
2D Toda lattice. This allowed to obtain a large class of explicit solutions to 
equation (15 .ip . the well-known hierarchy of multi-soliton solutions for exam- 
ple. The question of correct formulation of integrable boundary problem for 
2D Toda lattice (15 .11) was considered in [TJ] (see also [TB]). 

Here we are to describe a class of integrable boundary problems for (15 .ip 
and present examples of their explicit solutions. This procedure's ideology is 
similar to that applied above to KP equation (II .11) . so we will not get into 
details. 

The Zakharov-Shabat dressing method for (15 .ip consists of the following. 



uxy{X, Y, n) = w{X, Y,n-1)- w{X, Y, n), 



(5.1) 



\I^y(n) 



ux{n)'^{n) + -^{n + 1) = A°^(n), 
: -w{n - l)\&(n - 1) = L°\l/(ra), 



(5 .2) 



Introduce operator F so that 



oo 



(F$)(n) 



F{n,k)'^{k) 



(5 .3) 



fc=— oo 



[dx -A, F] = 0, [dy -L, F] = 0. 
Suppose that operator F admits triangular factorization: 



(5 .4) 



1 + F = {1 + Kr\l + K-), 



(5 .5) 
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where 



oo n— 1 

{K^)in) = Y,K{n,j)^{j), = J] ir-(n, j)^(j). (5.6) 

j=n j=-oo 

Then specify operators and of difference structure similar to that of 
operators L°, correspondingly by system 

idx-A')il + K) = il + K)idx-A'^) 

(dy - + ir) = (1 + ir)(9y - LO), • 

and obtain formula for new solution to 2D Toda lattice via initial solution 

u\n) =u%n) -log{l + K{n,n)) (5.8) 

The matrix K solves the difference analogue of Gelfand-Levitan-Marchenko 
(GLM) equation which arises from (15 .5p . (15 .6^ : 



K{n, m) + -F(n, m) + ^ K{n, k)F{k, m) = 0, m > n. (5 .9) 

k=n 

To launch the dressing procedure one needs to solve linear equations (15 .4p 
which look in the following way: 

dyFum = w{m)Fn,m+l - W{n - l)Fn-l,m- 

It is easy to see that system (15 .101) has solutions F^m = \l/(n)'Q/(m) where 
\l/(n) solves equations (15 .2p whereas \l/(n) solves the "conjugate" system 

^xin)=uxin)^in)-^in-l), 

§y(n) = w(n)§(n + 1). ^ ■ ^ 

In this case the desired solutions to GLM equation (15 .9p take the following 
form: 

K(n,m)^ "^'"'^'"l . (5.12) 

Now one has just to substitute (15 .12p into formula (15 .81) to obtain the 
"dressed" solutions, then in virtue of the Definition 1 select solutions which 
satisfy the integrable boundary constraint. 
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6 Integrable boundary problems for 2D Toda 
lattice. 

In systems of linear equations (15 .2p (15 .lip we replace variables (X, Y) — )■ 
{x,y) via functions X = X{x,y), Y = Y{x,y) and rewrite these systems in 
the following way: 

^.H = i {/22[-(/iiMx.(n) + fi2Uy{n))^{n) + ^{n + 1)] 
+f,,w{n~l)^{n-l)} 
= -i {/2i[-(/iiW. W + fi2Uy{n))^{n) + ^{n + 1)] 
+hMn-l)m{n-l)], 



= i {/22[-(/ii«.H + h2Uy{n))^{n) - §(n - 1)] 

^y(^) = -i {/2i[-(/iiW.(^) + h2Uy{n))^>{n) - §(n - 1)] 

+ 1)}, 



(6 .2) 



where fij are partial derivatives: /n = xx, /12 = /21 = xy, f22 = 
yy, A = /11/22 — /2i/i2- Suppose that solves the integrable boundary 
problem, therefore due to the Definition 1 the linear correlation 

= B{n)^{n) (6 .3) 

must hold on the boundary x = Xq or, in terms of variables X, F, on para- 
metric curve X = X{xo,y), Y = Y{xo,y). Substitute relation (16 .31) into the 
second equation of system (16 ■2p . using (16 .ip collect coefficients at functions 
\l/(n), ± 1) and obtain explicit form for multiplier B{n): 



(6 .4) 

x=xo 



and for the boundary constraint: 



= 0, 



ndy ^log^^ +Uy{n) + '^^{fnUx{n) + fi2Uy{n)) + dy log d{y) 

" (6 .5) 

where dy = d/dy, d{y) is arbitrary function. We impose also the following 
constraints: 

111 ^ 0, /21 ^ 0, A ^ 0. (6 .6) 
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7 Explicit solutions to boundary problems. 



Now we illustrate the method's effectiveness on a few examples. 
1. Take the following variable transformation: 

X{x,y) = e"(^')+''(^), Y{x,y) = e'=('^(^)-''(3')), (7 .1) 

where a{x), b{y) are arbitrary smooth functions, c is non-zero parameter. 
Then the boundary contour is presented by curve 



YX" = e2™(^'o) = D = const 



(7 .2) 



on (X, y)-plane. On this contour we have the appropriate boundary con- 
straint which ensues from formula (16 .51) : 



Ux{n) \x=xo = a'{xo) 



(l + c)n+- 



d'{y) 



b'{y)d{y)_ 

and the system of linear equations (16 .ip . (16 .2p reads as follows: 

^x{n) = a'{x)X^{n + 1) - i (w.(n) + ^{n) 
-ca'{x)Yu{n - - 1) 

%{n) = b'iy)X^{n + 1) - + ^^^^) ^{n) 

+cb'{y)Yu{n - - 1), 

^x{n) = ca'{x)Yu{n)^{n + 1) + i (n.,(n) + ^{n) 
-a'{x)X^{n - 1) 
= -cb\y)Yu{n)^{n + I) + + ^^^^) 

-b'{y)X^{n-l). 



(7 .3) 



(7 .4) 



(7 .5) 



Now we apply the integration procedure developed above to the boundary 
problem (15 .ip . (17 .2p . (17 .3p as d{y) = 0. Start from the following simple 
solution to problem fl5~T|) . (OD, (17^ : 

u\n) = {l + c)na{x). (7.6) 

Substitute (17 .6p into coefficients of systems of linear equations (17 .4p . (17 .5p 
and find their solutions: 
as c 7^ 1 



^(n) = exp (-f (1 + c){a{x) + 



2e 



(7 .7) 
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§H = exp (f (1 + c)(a(x) + b{y))) 



X Ell IjQj exp 



1-c V^i 



(7 i 



where kj, Ij, pj, qj, j = 1,2,...,N are arbitrary parameters. Now impose 
the integrabihty constraint 



^(n) = B(n)'^(n) 



\X=Xo ; 



where 

and obtain relations on parameters: 

1 . — h. . f.. — _Pip(T--c)a{xo) 

J J ' ' 

Then assume solution to GLM equation in the form (15 .12\i : 

"^(71)^(771) 



Kin, m) 



and substitute it into the formula 



u{n) = u^{n) — log(l + K{n, n)), 



(7 .9) 



(7 .10) 



(7.11) 



(7 .12) 



(7 .13) 



which yields the "dressed" solution to 2D Toda lattice that may be formally 
associated with A^-soliton solutions. In the simplest (A^ = 1) case we have 



K{n, n) 



V "exp 



^ )eS^*) 

c — 1 yp 



1 



1 - V 



-1 



which after substitution into (17 .13p gives the simplest "dressed" solution 
to the boundary problem. Now present the key formulae in terms of initial 
variables X, Y . 

Boundary constraint: 



Xux + Yuy = (1 + c)n 



(7 .14) 



boundary contour 



X'^Y = e2^«(^o) = D = const, 



(7 .15) 
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the simplest "dressed" solution to 2D Toda lattice is shown by formula 
dOSD, where 



= ^ (1 + c)n (^og X + ^ log F 



(7 .16) 



K{n, n) = 

^ - i.-" exp f-^ f-(l - +p(l--) X 



z/- 1 



c — 1 \p 



V 



, (7.17) 



•n i-c 
V = -—D — . 

c 



In the case c = 1 appropriate formulae take the following form: 
boundary constraint 



Xux + Yuy = 2n 



boundary contour 



XY = D, 

the simplest "dressed" solution to 2D Toda lattice is shown by formula 
(17 .13p . where 



u°(n) = n{logX + logF), 



(7 .18) 



K{n, n) 



-P 



2\-n 



Y 



1 

+ 



.XJ ' l+p2 

2. Now examine the variable transformation 



-1 



(7 .19) 



X = e'^^^) sin b{y), Y = ce'^(^) cos b{y) , 



(7 .20) 



which forms closed contour. Calculations which determine corresponding 
boundary problem and its exphcit solutions are similar to those mentioned 
above, so we present here just final outcome in the initial variables: 
boundary constraint 



/Y X " 
Yux + c^Xuy + ^ ^ " °' 



boundary contour 
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y2 

the "dressed" solution is presented by formulae 

= -nlog(Xy), 

K{n, n) = 



+ 



1 + 5 



(7 .21) 



-1 in the 



p is arbitrary parameter. 

3. At last we have the example of regular solution. Put c 
formulae (17 -ip - dT .17p and obtain: 

boundary constraint 



Xux — Yuy = 0, 

boundary contour 

Y/X = e2"(^°) = D. (7 .22) 

The initial solution is chosen as a constant: u^{n) = vP and 



K{n, n) 



z/- 1 



V ""exp 



(7 .23) 



where v = j^. The formula fl7 .13p yields solution regular as < z/ < 1. 
This solutions presents an ordinary soliton that travels along the straight 
boundary line (17 .221) . so no reflection effects occur. 



8 Conclusion. 

In this paper we show that: 

1) integrable in force of the Definition 1 boundary problems for Kadomtsev- 
Petviashvili equation and 2D Toda lattice can be solved explicitly. We develop 
the method and show its effectiveness on concrete examples 
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2) the set of integrable boundary problems that consist of boundary coun- 
tours and appropriate boundary constraints is large enough. This fact is il- 
lustrated by formulae dOJ, dOD, §IM- 

The author thanks Professor I.T. Habibullin for valuable consultations. 
The investigation was partially supported by grants of the Program of RAS 
Presidium "Fundamental problems of nonlinear dynamics" and Russian Foun- 
dation for Basic Research 10-01-00088-a. 



15 



References 

[1] B.B. Kadomtscv and V.I. Petviashvili// Soviet Math. Dokl., v. 192, n. 

4, p. 753-75 (1970) 

[2] Mikhailov A.V.// JETP Letters (1979), v.39, n.7, P. 443-448 

[3] S.P. Novikov, S.V. Manakov, L.P.Pitaevski and V.E. Zakharov. Soliton 
theory. Nauka, Moscow 1980 

[4] V.S. Dryuma// JETP Letters (1974), v. 19, n. 12, n. 753-755 

[5] V.E. Zakharov and A.B. Shabat// Funkz. anahz i ego pril. (1974), v. 8, 
n. 3, p. 43-53 

[6] E.K. Sklyanin// Funkz. anahz i ego pril 21:2 (1987), p. 86-87 

[7] I.T. HabibuUin// Theor. Math. Physics, 91:3 (1992), 363-376 

[8] I.T. HabibuUin and E.V. Gudkova// Funkz. analiz i ego pril. 38:2 (2004), 
71-83 

[9] I.T. HabibuUin and E.V. Gudkova// Theor. Math. Physics, 140:2 (2004), 
230-240 

[10] V. Adlcr, B. Gurel, M. Gurses and I.T. HabibuUin //J. Phys. A: Math. 
Gen. 30 (1997), p. 3505-3513 

[11] Mikhailov A.V.// Physica D (1981), V.3, n.1,2, P.73-117 

[12] K. Ueno and K. Takasaki, Toda Lattice Hierarchy, Adv. Stu. in Pure 
Math. 4 (1984) 1 

[13] Nakamura A. // J. of the Phys. Soc. Japan, (1983), V.52, n.2, P.380-387 

[14] M.V. Babich, V.B. Matveev and M.A. Sail// Zapiski nauch. seminar. 
LOMI (1985), V.145, p. 34-45 

[15] Gurses M., HabibuUin I.T. and Zheltukhin K.// J. of Math. Physics 
(2007), V.48, 102702 

[16] HabibuUin I.T. and A.N.Vil'danov, Integrable Boundary Conditions for 
Nonlinear Lattices, CRM, Proceedings and Lecture Notes, Amer. Math. 
Soc, Providence, v. 25, 2000, p. 173 - 180. 



16 



